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A variety of important physical problems cannot be addressed within the framework of fixed-order perturbation 
theory, the most widely used calculational scheme in the continuum. This is often the case within Quantum Chro- 
modynamics (QCD), when large disparities in the physical scales involved result in a complicated interplay between 
perturbative and non-perturbative effects. Similar limitations appear when physical kinematic singularities, such as 
resonances, render the perturbative expansion divergent at any finite order, or when perturbatively exact symmetries 
prohibit the appearance of certain phenomena, such as chiral symmetry breaking or gluon mass generation. In such 
cases one often resorts to various reorganizations of the perturbative expansion inspired from scalar field theories or 
Quantum Electrodynamics (QED), supplemented by a number of auxiliary physical principles. When studying the 
interface between perturbative and non-perturbative QCD for example, one finds it advantageous to use concepts 
t-H 1 familiar from QED, such as the effective charge, in conjunction with dispersive techniques and analyticity properties 
of the S'-matrix [|lj. In addition, in the field of renormalon calculus, one studies the onset of non-perturbative effects 
from the behaviour near the QCD mass-scale of judiciously selected infinite sub-sets of the perturbative series [||. 
Similarly, the extension of the Breit-Wigner formalism to the electro- weak sector of the Standard Model necessitates 
' a non-trivial rearrangement of the perturbative expansion Q ; an analogous task must be undertaken when studying 
t— ( [ various aspects of finite temperature QCD |I| , as well as mass generation, both in 3-dimensional field-theories || and 
ON in QCD as a prelude to a systematic truncation scheme for the Schwinger-Dyson series. 

One of the main difficulties encountered when dealing with the problems mentioned above is the fact that several 
physical properties, which are automatically preserved in fixed-order perturbative calculations by virtue of powerful 
Oh, field-theoretical principles, may be easily compromised when rearrangements of the perturbative series, such as re- 
summations, are carried out. These complications may in turn be traced down to the fact that in non-Abelian gauge 
theories individual off-shell Green's functions (rt-point functions) are in general unphysical. 
■ It turns out that this last problem can be circumvented by resorting to the method known as the pinch technique 
! (PT) dJ^]. The PT reorganises systematically a given physical amplitude into sub-amplitudes, which have the 
same kinematic properties as conventional ro-point functions, (propagators, vertices, boxes) but, in addition, are 
endowed with desirable physical properties. Most importantly, at one- loop order (i) are independent of the gauge- fixing 



, — i 

X 



■ parameter; (ii) satisfy naive, (ghost-free) tree-level Ward identities, instead of the usual Slavnov- Taylor identities, (iii) 
contain only physical thresholds and satisfy very special unitarity relations |^|| (iv) coincide with the conventional 
n-point functions when the latter are computed in the background field method Feynman gauge (BFMFG) jnj. These 
properties are realized diagrammatically by exploiting the elementary Ward identities of the theory in order to enforce 
crucial cancellations , and make manifest intrinsic properties of the S-matrix, which are usually concealed by the 
quantization procedure. 

The important question which arises is whether the PT algorithm may be extended beyond one-loop, leading to the 
systematic replication of the aforementioned special properties of the PT effective n-point functions to higher orders 
jjJI . In this Letter we will show that the PT can be generalized to two loops by resorting exactly to the same physical 
and field-theoretical principles as at one-loop. 

We start by briefly reviewing the one-loop case. Consider the S-matrix element for the quark (w)-antiquark (u) 
scattering process u{P)u(P') — > u(Q)u(Q') in QCD; we set q = P' — P = Q 1 — Q, and s — q 2 is the square of the 
momentum transfer. It is convenient to work in the renormalizable Feynman gauge (RFC); this constitutes no loss 
of generality, since the full S'-matrix is independent of the gauge-fixing parameter and gauge-fixing scheme. One first 
decomposes the elementary tree-gluon vertex Ta}iu{q,Pi,P2) as follows 0: 



^alAi,v\,P2) = [{pi -P2) a g^v + lq v g ai i - ^9 av \ + [P2vg af . 
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This decomposition assigns a special role to the q-leg, and allows to satisfy the Ward identity 

<fTf a ^piiP2) = (pl-pl)g^ (2) 

where the right-hand-side is the difference of two- inverse propagators in the Feynman gauge, and vanishes on shell, i.e. 
Pi = P2 = 0- Notice that the first term in Tp v is a convective vertex, whereas the other two terms originate from 
gluon spin or magnetic moment. ^p\ lnj {q_,Pi,P2) coincides with the BFMFG bare vertex involving one background 



(q) and two quantum (^1,^2) gluons [ji0|. 

We then carry out the above decomposition on the three-gluon vertex appearing inside the non-Abelian graph 
contributing to the one-loop quark-gluon vertex . The result of this is two-fold: First, the action of the longitudinal 
momenta p^ = — fc M , p^ = {k — q) u on the bare quark-gluon vertices and respectively, triggers the elementary 
Ward identity of the form ft = ( ft+ fQ — m) — (fQ — m). The first term gives rise to the pinch contribution Vp^ a (q) 
given by Vp^ a {q) — 2q 2 Ca J[dk][k 2 (k + q) 2 ]^ 1 g aa , where g is the gauge coupling, Ca is the Casimir eigenvalue of the 
adjoint representation, and [dk] = p 2t ^J^ d , with fi the 't Hooft mass; the second term vanishes on-shell. Second, the 
part of the graph containing together with its Abelian-like counterpart defines the PT one-loop quark-gluon 

vertex f^(Q,Q'), which satisfies the QED-like Ward identity q a T^\Q,Q') = - £ (1) (<2')> where £W is 

the PT one-loop quark self-energy. The propagator-like parts extracted from the vertex are cast into the form of a 
genuine self-energy by setting Ilp^(g) = Vp 1 ^ (q) a <jtp(q), where t^{q) — q 2 g l _ LV ^q^q v ', thus, the resulting one-loop PT 
self-energy reads II^ 1 ! (q) — II^ 1 ] (q) + Hp^p (<?) ■ Carrying out the one- loop integrations one finds || that the prefactor 

in front of the logarithm of TL^l(q) is (11/3)CU, i.e. the coefficient of the one-loop j3 function for quark-less QCD. 

For the two-loop case, one considers the two-loop S-matrix element for the aforementioned process uu — ► uu in the 
RFG, and focusses on the two-loop quark-gluon vertex T y a '(Q, Q'). The Feynman graphs contributing to r„ (Q, Q ) 
can be classified into two sets, (a) those containing an "external" three-gluon vertex i.e. a three-gluon vertex where 
the momentum q is incoming (Fig.l). (b) those which do not have an "external" three-gluon vertex. This latter set 
contains either graphs with no three gluon vertices (abelian-like) , or graphs with three-gluon vertices whose all three 
legs are irrigated by virtual momenta, i.e. q never enters alone into any of the legs. Carrying out the decomposition of 
Eq. ([!]) to the external three-gluon vertex of all graphs belonging to set (a) , leaving all their other vertices unchanged 
, the following situation emerges: 



ri 2 >(Q,Q') = lf)(Q,Q') + ^(q)^ > + ~n^( 9 )(^)r«(Q, Q') , (3) 



with 



V ( pIM = -h 



(l - q) p + (2J 2 + h) 9o 



k, -k - i)T^ x »{l, k, —k — tj — 2k a (k + tU , (4) 



where h = I (k + iy 2 {k + I - q)- 2 , I 2 = I Q (k + q)- 2 , I 3 = I (k + I)- 2 , h = I £- 2 (k + i)~ 2 , with il 
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g 4 C\[£ 2 (£ — q) 2 k 2 ] 1 , and the two-loop integration prefactor (p 2e ) 2 j J 7§^j§~tJ has been suppressed. Ya (Q,Q') 



(2) 

is the two- loop BFMFG quark-gluon vertex, (g) the propagator- like part, and the third term on the right-hand 
side is the necessary contribution for converting the one-particle reducible part of the two-loop S'-matrix element 
rL° ) (^)n^(q)(^)rJ ; 1) (g,Q') into rL° ) (^)n2( 9 )(^)f^ 1) (Q,Q')- Eq.® is a non-trivial result, since there is no 
a-priori reason why the implementation of the decomposition of Eq. (|l|) should only give rise to terms which can be 
interpreted in the way described above. In fact, individual diagrams, or even natural sub-sets of diagrams such as the 
one-loop three-gluon vertex nested inside the two-loop quark-gluon vertex, give in general rise to contributions which 
do not belong to any of the terms on the right-hand side of Eq. (0) . It is only after all terms have been considered that 
the aforementioned crucial cancellations become possible. Finally, the counterterms of r„ (Q,Q') must be correctly 
accounted for Q. fi 2) (Q, Q 1 ) satisfies the QED-like Ward identity q a f L 2) (Q, Q') = £ (2) (Q) - £ (2) (<3'), wher e £ (2) is 
the two-loop PT quark-self-energy. X^ 2 ) is identical to the conventional £( 2 ) in the RFG (and the BFMFG), exactly 
as happens at one-loop. 

To construct the two-loop PT gluon self-energy n^,i(g), one must append to the conventional two-loop self-energy 
ngfa) the term lif^(q) = v£,(«)tg(g) together with the term iR^(q) = 10^%) + {q)n^{q) 
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originating from converting a string of two conventional one-loop self-energies into a string of two one-loop PT self- 
energies H . One can show by means of a diagram-by-diagram mapping that the resulting H^l (q) is exactly identical 
to the corresponding two- loop self-energy of the BFMFG, and that this correspondence persists after renormalization 
jl4]]. Notice that the presence of the term Rp a g(q) is crucial for the entire construction, and constitutes a non-trivial 
consistency check of the resummation mechanism first proposed in ||. An immediate consequence of the above 

r(l)/ 



correspondence is that the coefficient in front of the leading logarithm of H a p(q) is precisely the coefficient of the 
two-loop quark-less QCD j3 function [Q, namely (34/3)C^ |Tsj] . As a result, one may extend to two-loops the one- 
loop construction of a renormalization-group- invariant effective charge presented in [ |l6| , leading to the unambiguous 
identification of the conformally-(in)variant subsets of QCD graphs Finally we note that, exactly as happens at 
one-loop, the two-loop PT box-graphs are simply the conventional ones in the RFG (and are equal to the ones in the 
BFMFG). 

As has been explained in detail in the one- loop PT n-point functions satisfy the optical theorem individ- 

ually. To verify that one starts with the tree-level process u(P)u(P') — > g(p\) + g{p2), whose S-matrix element 
we denote by 7^; then, one considers the quantity T^P^ (pi)P uv [p^T^yi, where P M „(p, rj) = -g^ v + [rj^Pv + 
VvPn) I TP + rfp^Pv I (vp) 2 i with rj an arbitrary four-vector. One proceeds by first eliminating the Ip^ (g,pi,p2) part 

of Ta}iv{q,Pi,P2), which vanishes when contracted with the term P w ' (pi)P u "' (pz). Then, the longitudinal parts of 
the P MAl '(pi) and P VV '{P2) trigger a fundamental cancellation involving the s- and t- channel graphs, which is a 



consequence of the underlying BRS symmetry Jl5| . Specifically, the action of pi M on the Ij?^ gives 

PlTFauui&Pl'Pz) = t a v{q) + (pi -P\)g a v + (P2 -P\)a.P2v \ 



(5) 



the first term on the right-hand side cancels against an analogous contribution from the i-channel graph, whereas the 
second term vanishes for on-shell gluons. Finally, the term proportional to pi v is such that (i) all dependence on i] 
vanishes, and (ii) a residual contribution emerges, which must be added to the parts stemming from the g^Qw' part 
of the calculation. Then one simply defines self-energy/ vertex/box- like sub-amplitudes according to the dependence 
on s = (pi + P2) 2 and t = (P — pi) 2 , as in a scalar theory, or QED. The emerging structures correspond to the 
imaginary parts of the one-loop PT effective Green's functions, as one can readily verify by employing the Cutkosky 
rules; in fact the residual pieces mentioned at step (ii) above correspond precisely to the Cutkosky cuts of the one- loop 
ghost diagrams. The one-loop PT structures may be reconstructed directly from this tree-level calculation, without 
resorting to an intermediate diagrammatic interpretation, by means of appropriately subtracted dispersion relations. 

The same procedure must be followed at two-loops; the only difference is that one must now combine contributions 
from both the one-loop 5-matrix element for the process u(P)u(P') — > g{pi) + g{p2) and the tree-level S'-matrix 
element for the process u(P)u(P') — > g(pi) + g{P2) + 9(^3) ■ The non-trivial point is that the one-loop S'-matrix 
element must be cast into its PT form (as shown in Fig 2a.) before any further manipulations take place. Notice 
that the same procedure which leads to the appearance of 11(g) leads also to the conversion of the conventional 



one-loop three-gluon vertex ^ap,v(q,Pi,P2) into Tp* a (q,pi,p2), which is the BFMFG one-loop three-gluon vertex 

with one background (q) and two quantum (pi, p^) 
the following Ward identity 



ft is straightforward to show that Tp a „ v (q,Pi,P2) satisfies 



rt&^fopx.pa) = nWfa) - n«(p 2 ), 



(6) 



which is the exact one- loop analogue of the tree-level Ward identity of Eq (g); indeed the right-hand side is the 
difference of two one-loop self-energies computed in the RFG. In order to extend to the next order the dispersive 
construction outlined above, one needs the following Ward identity 



p^fL, = *n«( Q ) - fflWfe) + A«ts(«) 



(7) 



with 



(k - pi)'4°)> 2 , k, -k - V2 ) - (k+ P2 )uK 



P2k p Tpl ap {q,k+p2, -k + pi)-p 2 - (k -pi)(2k +P2 -Pi) a 



2B(q)+B(p 1 ) 



gua 
I 



B(pi) + B(p 2 ) 



(8) 
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J = \ig 2 C A [k 2 {k - P i) 2 (k +P2) 2 ]- 1 and B{p) = g 2 C A j [dk][k 2 {k + q) 2 }- 1 . Eq. (@) is the one-loop analogue of Eq. 
(||) . The one-loop version of the fundamental BRS-driven cancellation will then be implemented; for instance, the 
first term on the right-hand side of Eq. (]?]) will cancel against analogous contributions from the graph of Fig. 2a2, 
whereas all remaining terms proportional to t aa (q) will cancel against contributions from the i-channel graphs of Fig. 
2a 3 

The same construction must then be repeated for the tree-level process uu — * ggg, whose tree-level S'-matrix element 
we denote by T^vp ; again, the s-channel graphs (Fig. 2b) must be rewritten in such a way that when contracted with 
q only terms proportional to p 2 emerge, but no transverse pieces, exactly as in Eq.(Q). This is accomplished by simply 
carrying out the decomposition of Eq.(Q) only to the vertices where q is entering; then the contributions originating 
from the Tp v parts eventually vanish when contracted with the polarization tensors P w (p\)P vv {p2)P pp ips)- 
Acting with the longitudinal parts of the polarization tensors on the T^vp one must first carry out the corresponding 
BRS s — t channel cancellation, and pick up automatically the correct ghost parts. Notice in particular that this 
procedure gives rise to the ghost structure given in Fig. 3c of jL7|, which has only three-particle Cutkosky cuts, and 
does not exist in the conventional formulation. 

Adding the s-channel terms together the total propagator-like part emerges; it is proportional to (34/3)C^g 2 , as it 
should. Notice that the result is infrared finite, by virtue of crucial cancellations between the the one-loop uu — > gg 
and the tree-level uu — > ggg cross-sections. The most direct way to verify that is by exploiting the one-to-one 
correspondence between the terms thusly generated and the Cutkosky cuts of the BFMFG two-loop self-energy; the 
latter are infrared finite since they effectively originate from a single logarithm. 

In conclusion, we have shown that the same physical principles, and, evidently, the same procedure used at one-loop, 
lead to the generalization of the PT to two- loops. In particular , the known correspondence between PT and BFMFG 
persists. It would be interesting to explore its origin further, and establish a formal, non-diagrammatic understanding 
of the PT. 
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FIG2. The one-loop s-channel corrections to uu — » gg, and the tree-level s-channel uu — » ggg 
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